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In this paper, colorless bilocal fields are employed to study the large N limit of 
both fermionic and bosonic vector models. The Jacobian associated with the change of 
variables from the original fields to the bilocals is computed exactly, thereby providing 
an exact effective action. This effective action is shown to reproduce the familiar per- 
turbative expansion for the two and four point functions. In particular, in the case of 
fermionic vector models, the effective action correctly accounts for the Fermi statistics. 
The theory is also studied non-perturbatively. The stationary points of the effective ac- 
tion are shown to provide the usual large N gap equations. The homogeneous equation 
associated with the quadratic (in the bilocals) action is simply the two particle Bethe 
■ Salpeter equation. Finally, the leading correction in -k is shown to be in agreement 

with the exact S matrix of the model. 
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I. INTRODUCTION 

> 

■ One of the most important problems facing high energy physics today is the construction of suitable approximation 
techniques that will allow an analytic understanding of the long distance properties of non Abelian gauge theories. 

■ The large N expansion remains one of the most promising techniques proposed to date. The leading order of this 
approximation, the master field, is given by the sum over planar diagrams in matrix models [Q. It has up to now been 
impossible to compute this sum in closed form, except for simple systems §]. For this reason, not even the leading 

I term in the expansion, in four dimensions, has been computed. In principle, if the leading term could be obtained, 
""^ ■ one would set up a perturbation theory around asymptotic states consisting of "mesons" and "hadrons" - colorless 
bound states of the quanta of the matter fields appearing in the Lagrangian. In other words, fluctuations about this 
^ . master field yield the mass spectrum of the theory p . 

\ Vector models by contrast, have soluble large N limits. Indeed, the large N limit of a number of vector models has 
^2 ' been studied using an auxiliary field ^ , [|J , which is not (classically) dynamical. Formally eliminating the auxiliary 
field, the original model under consideration is regained. The advantage of the auxiliary field is that it is chosen 
so that it contains no implicit N dependance. Thus, integration over the original field variables yields an effective 
theory in which all N dependence is fully explicit. There are however, at least two serious objections to this approach. 



Firstly, it is only viable for Lagrangians consisting of quadratic plus quartic terms. Secondly, the local auxiliary field 
" ' is a poor substitute for the composite (bilocal) meson field. 

Jevicki and Levine || have constructed the master fields for a large number of bosonic vector models, working 
directly with equal time bilocal meson fields, thereby overcoming the difficulties mentioned above. They rewrite the 
Hamiltonian in terms of these bilocal fields, obtaining an effective theory of mesons in which all N dependance is 
fully explicit. This change of variables induces a non trivial Jacobian [Q, which is obtained by imposing hermiticity 
on a suitably scaled Hamiltonian. This approach is not easily implemented for fermionic vector models: it is a well 
known fact that the Dirac hamiltonian describes constrained dynamics. These constraints obscure the Hamiltonian 
approach. For this reason, the large N limit of fermionic vector models was originally studied using a pseudospin 
formalism ||. Although this equal time approach does provide a possible starting point for a systematic expansion 
in 1 /N, the presence of the pseudospin constraints obscures the relationship between the pseudospin variables and 
quantities of interest. The possibility of reexpressing a fermionic vector theory in terms of unequal time bilocals has 
been recently considered by Cavichi et al. || who obtained the leading large N form of the effective action directly in 
terms of these bilocals, and were therefore able to reproduce t'Hooft's QCD2 equation Jl(| for the meson spectrum. 

In matrix models, higher order corrections in have been systematically computed in the case of the c=l matrix 
model string theory |pd[| directly in terms of invariant variables. These corrections have been shown to be in complete 
agreement with the exact solution of the model. In this case the effective propagator is simply that of a massless scalar 
and it also turns out that the effective Hamiltonian is cubic. The invariant variables of the matrix model correspond 
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to the bilocals in vector models, so the question that we address in this article is whether systematic higher order 
corrections in can be computed directly in terms of the bilocal variables for vector models. 

Higher order calculations have been performed for a Gross Neveu model with scalar-scalar interactions, using the 
auxiliary field approach fll2| , or directly in terms of the original fermionic fields Jl3[ as a check of the expansion, 
given that the exact S matrix of the model is known |l4j]. 

In this paper, we develop the necessary formalism needed to systematically compute higher order corrections in j*, 
by obtaining an exact effective action in terms of unequal time bilocal fields. These bilocal fields are time ordered 
and we use them to study both bosonic and fermionic vector models. 

We obtain the exact effective action by explicitly calculating the full Jacobian associated with the change of variables 
from the original fields to the time ordered bilocals. This is achieved by requiring that the Schwingcr Dyson equations 
derived directly in terms of the bilocal variables agree with the equations derived in terms of the original fields [ p"5[ . 
This provides a unified approach for both bosonic and fermionic models. The bosonic Jacobian agrees with the result 
obtained from collective field theory 0. The result for the fermionic Jacobian is new. It is remarkable that within 
the context of a path integral quantization such a Jacobian exists and is exactly computable. 

This effective action is nonlocal in time but in the context of the functional integral quantization it generates all 
the "colorless" correlators of the quantum theory. For fermionic theories, since the effective action is exact and also 
valid for N = 1, it provides an exact rewriting of the original theory (although for N = 1 of course there is no 
small expansion parameter). We therefore expect this approach to have interesting applications in condensed matter 
type fermionic systems. The formalism is also likely to be particularly relevant for problems in which the spectrum 
is expected to consist only of singlets under the global invariance of the theory, since it then directly describes the 
dynamics of the physical degrees of freedom. Finally we wish to emphasize that this formalism provides a (nonlocal) 
bosonization valid for arbitrary dimensions. Although many of these aspects are now under study the main purpose 
of this article is to obtain the effective action and confirm its validity perturbatively. 

So, perturbatively in both -h and the coupling constant, we consider the linear sigma model and the Gross Neveu 
model (with scalar-scalar interactions) and verify that typical correlators are in agreement with their Feynman ex- 
pansions, and in particular verify that the statistics is properly taken into account. This provides a strict test of the 
validity of our effective action. We go on to describe how to develop a systematic expansion of the effective action 
and explicitly establish as the perturbative parameter of the expansion. 

The paper is organized as follows: in section II the effective action, connection to the collective field theory and 
check of the perturbative expansion in both and the coupling constant is carried out. In section III we do similarily 
for fermionic vector models. In section IV we make contact with the (nonperturbative) analysis of the spectrum for 
Gross Neveu and Nambu-Jona-Lasinio type model by obtaining the propagator of our effective field theory. We show 
further that the homogeneous equation for the quadratic fluctuations of our effective field theory is a Bethe-Salpeter 
equation for these models. In section V we obtain all diagrams contributing to the correction to two particle 
scattering from the effective field theory and show that they are in precise agreement with the diagrams considered 
in previous comparisons of the j? expansion with the exact S matrix p3| . 

II. BOSONIC VECTOR MODELS AND THE RELATIONSHIP TO THE COLLECTIVE FIELD THEORY. 



In this section we will consider U(N) invariant bosonic vector models. By this we mean any theory of N complex 
scalar fields 4> a {x), a = 1, 2, ...N (x is a d dimensional parameter) with an action S that is invariant under the global 
symmetry cj) a — > <f)' a = U ab (j) b , <ft* a — > cj)'* a — U ab (j)* b with U ab an arbitrary U(N) element. Furthermore we assume 
that all coupling constants of the theory have been rescaled appropriately so as to yield a systematic 1/iV expansion. 
This is equivalent to the statement that a rescaling of the fields exists under which S — > NS. A specific example that 
will be examined in detail later in this section is given by the action 



although our discussion applies to arbitrary invariant actions S as described above. These include (non local) effective 
actions resulting from the explicit integration of intermediate degrees of freedom such as QCDi. This will be the 
object of further study in another communication Jig], 



A. Schwinger Dyson equations and Jacobian. 




(1) 
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We will assume further that one is only interested in time ordered product expectation values of invariant operators 
(generically referred to as invariant correlators) 

< fM > s <n^)>= ^-«;<-> 

4 — 1 

with 

a{x,y)=4>* a {x)<f> a {y) (3) 

It is straightforward to obtain a set of Schwinger-Dyson equations for the invariant correlators. These follow from the 
identity 



= / D(j>*D(f) 



5(j) a {x) 



<t> a (y)F[ay b (4) 



yielding 



< N5 d {x - y)F[a] > + < > < W^L^M >= ( 5 ) 



The above set of equations involve only invariant correlators as it will be shown in the following. 

The philosophy that we wish to adopt here is that there is a Jacobian J associated with the change of variables 
from the original variables 4> a to the invariant variables (||) inside the functional integral that correctly yields all time 
ordered product expectation values of invariant operators [Q. In other words, 

J D(j)* Dcj)F[o]e lS = J DaJF[a]e lS (6) 

We now follow the procedure described in reference [jl5| to obtain a differential equation for the Jacobian from the 
loop equations This results from the identity: 

= I Da [ d d z—^-{a{z,y)JF[a]e lS ) (7) 
J J 6a(z,x) 



which implies^ 

5a{z,x) J ' 5a(z,x) J ' 5<j(z,x) 



<S d (x-y)S d {0)L d F[a} + I d d za(z,y)-^—F[a] + f d d za(z,y) +i I d d za{z,y) 6S - F[a] >= 0. 



(8) 



We can now use the chain rule 



5 f d d z f d d y S ^l4_±_ = ( d d z ^ {z) _±_ (9) 



S<j) a (x) J J S<j) a (x) Sa(z,y) J Sa(z,x) 

in equation (j^) which is then equivalently written as 

f S F f S S 

< NS d (x - y)F[a] > + < J d d za(z, y)j^^ > +i < J d d za(z, y) F[a] > = 0. (10) 

Requiring agreement of this last equation with (^|) for arbitrary F[a] we obtain 

d d za(z, v)t^t = (N- L d 5 d (Q))5 d (x - y). (11) 
oa(z, x ) 



x We use L d to denote the volume of the system we are studying, i.e. L d = J d d x. 
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J is independent of the action, as it should be. The solution to this equation is: 

InJ =(N - L d 5 d (0))Trlna (12) 
where the trace is in functional space. For models with "flavor" degrees of freedom (\> a a , a = 1, 2. ..to the invariants are 

a a0 (x,y) = ^(^(y) (13) 

It is straight forward to generalize the above analysis and to show that the Jacobian is given by 

InJ = (N - mL d S d (0))trTrlna (14) 

where the trace is now taken in both functional and flavor space. This result generalizes for an arbitrary number of 
dimensions the result obtained in (ItJ for d — bosonic vector models. 

B. Connection with the collective field theory 

The point of view that the large N limit can be understood in terms of a change of variables to invariant variables 
or subspaces was successfully exploited by Jevicki and Sakita fij for a large class of models. In the first of @ these 
authors are able to obtain the form of the effective Hamiltonian acting on the reduced invariant subspace in terms of 
the Jacobian associated with the new inner product measure. Remarkably the equation satisfied by the Jacobian 
can be obtained by the simple requirement that the effective Hamiltonian must be explicitly ^] Hermitian and is given 
by 

Vo fr n 6lnJ - U (c\ y gXgigO (15) 

\ n{c > c) sW)- " {c) - \ "WT (15) 

In this equation C and C index the invariant variables. In our case 

V(x,y; X ',v') = J d d -^ I ^ ) a*( X ,y)-^a(x',y')=6' 1 -\y-y')a(x', X ) (16) 

and 

»> = I ^ v° (*)V(*) a{x > v) = md ~ 1(x - y) - (i7) 

The above equations only involve equal time correlators as it is appropriate for a Hamiltonian approach. 

Jevicki and Sakita, in the second of j?J have shown that exactly the same equation must be satisfied in a functional 
integral description provided the invariants are appropriately labeled i.e. in our case if x and y are d dimensional 
spacetime points. With this proviso, it is straightforward to show that once equations (|l6|) and ([Tt]) are used in 
equation (fL5[) we reproduce the differential equation (jll|) for the Jacobian derived in the previous subsection. 



C. Perturbative Check 

The full Jacobian @ is not new; it has certainly been written down in ||. However in all applications that we are 
aware of only the leading large N term in the Jacobian has been used to obtain the leading contribution to the free 
energy, leading time ordered correlators and spectrum It is actually not clear in what sense the second term in 
the Jacobian is " subleading" due to infinities appearing in it and it is not at all obvious how that it can help generate 
systematic 1/.ZV corrections. It is the purpose of this subsection to show that this is indeed the case, perturbatively 
in the coupling constant. 



2 With respect to the trivial measure. 
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1. Exact Effective Action 



Any invariant correlator can be calculated from 

f DaJe lS F\a] f Dae lS -ff F\a] , x 

< F\a] >= J— r = J — f — 18 

where 

S eff = -ilnJ +S = -iNTrlna + S + iL d 8 d (0)Trlna (19) 

In order to exhibit explicitly the N dependence we rescale a — > No under which, as explained in subsection A, 
S — > NS. Throughout the rest of this section this rescaling is always implicit. Therefore we have: 

S eff = -iNTrlna + NS + iL d S d {0)Trlna = NS + St. (20) 

We see that as N — > oo the leading (" classical" ) configuration is determined by 

^Uo=0 (21) 



We can now perturb about a by letting 



a(x, y) = a°{x, y) + -^=Mx, y) (22) 
V N 



and expanding S e ff as 

S eff = NS (cr°) + S 1 (a°) + -B 2 + -A 2 + £ -=[ - l ^J^B n+2 + iL d 6 d (0) [ 



B r , + —A 



n+2\ 



where 



and|^ 



W L (n + 2) ^ v ' n n (n + 2)! 

(23) 



A n = I dxi... I dx n I dyi... I dy n — / S . A c oT]{xi, y x ) r](x n ,y n ) (24) 

oa{Xi,yi) da(x n ,y n ) 



B n = Tr(a° l v ) n . (25) 



We obtained an effective action with an infinite number of vertices as it is to be expected from any loop expansion. 
It should however be remembered that in order to calculate any diagram to a given order of -7= only a finite number of 
vertices need to be included. Notice that the subleading term of the Jacobian (|l^) induces tadpole type interactions. 
In the case of c = 1 strings where the effective Hamiltonian consists of a cubic and a tadpole interaction this tadpole 
interaction has been shown to be essential for an agreement with exact results |Oj. For d = vector models similar 
arguments have been presented in reference |l7|]. We will see in the following that the tadpole interaction will be 
essential to obtain agreement with a perturbative Feynman analysis of the (cf>* a (f> a ) 2 theory. 



3 The expansion of Trlna in this fashion is justified by the fact that translational invariance requires <r° to be diagonal in 
momentum space. 
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2. Perturbative Results for the {(f>* a <j> a ) 2 theory. 
We will now be considering in detail the theory defined by the action (TO) 



S = I d d x(d^* a d^ a - m 2 r a <p a ~ > Q ) 2 ) (26) 



We will need the following perturbative result for the two point function: 



• ""i.rk<"(/,) = / ^G 2 {jp) (27, 



where 



G 2 (p)=N( -2*£ Q +4(^) 2 Q Q + 4 (^) 2 Q +0(g 6 )) 

+(-2# Q +srf) 2 Q Q + 8 (^) 2 8 +4(^) 2 ^^+0(. g 6)) +0( ^) 



8 ; 1 /; 1 ^\ N' 

(28) 



For the connected piece of the four point function 



where 



Gr n (Pi,P2,P3,P4)=iV(^ /\-fi -S 7\ +0(.g 6 )) + 0(l). 

(30) 

In the above expressions the diagrams have the standard Feynman interpretation (without symmetry factors and 
coupling constants). For instance 



if k i i 



(31) 



We assume that all diagrams have been suitably regularized. Our aim in this section is not to discuss issues 
related to renormalization, but to confirm that our effective theory systematically reproduces the l/N expansion, 
perturbatively in g 2 . 

3. The Leading Order. 

With a translationally invariant ansatz 

<r(x,v) = J -^e^^a(p) (32) 
It is straightforward to verify that the solution to equation (QW) a (p) satisfies 



*°(P) = n rj^ — (33) 



6 



This is the familiar gap equation. Iterating this equation to 0{g A ) one obtains 



<7°(p) 



^jL+^fmi +4(-g^) 2 JL +0(g 6 ) 



(34) 



in agreement with the O(N) term in equation (p8|). (We recall that in our effective theory the fields are rescaled by 
factors of N.) 



4- Effective Field Theory Spectrum 



The leading order <t° can be used, at the level of quadratic fluctuations, to obtain the mass spectrum of the 
theory Perturbing about the leading order, as has been described, we obtain the (leading) quadratic action 

i 1 

1 f 2 

= - dx 1 dx 2 dx 3 dx 4: [-^-S(x 1 - x 2 )S(x 1 - x 3 )5(x 1 - x 4 ) + id (x 1 ,x 2 )a° (x 3 ,x 4 )]i](x 4: ,x 1 )ri(x 2 ,x 3 ). (35) 



After introducing a symmetric Fourier transform 

d d P i f d d P2 



(2n) d / 2 J (27r) d / 2 
we can write the quadratic action in momentum space 



(36) 



S * = ~2 J dpidp2dpadpi l j(2m¥ S ( pl ~ P2 +P3 + iS ( pi ~ Pz) S (Pz ~ P^v (Pi) a ° (P3)]v(P±,Pi)v(P2,P3) (37) 

This action determines the propagator A(pi,p 2 ,p 3 ,p 4 ) of the a field. Conventional arguments show that 
A(pi,P2,P3,Pa) satisfies the equation: 



dpidp 2 



g 2 S(ki +pi - k 2 -p 2 ) 



4 (2ir) d 
X-A(Pi,P2,P3,P4) = io~(k 2 -Pi)S(ki ~p 3 ) 



+ iS{h - p 2 )<r° (pi)tf(pi - k 2 )a° (p 2 ) 



(38) 



Now, the physical interpretation of A(pi,p 2 ,p 3 ,p 4 ) is as follows: this propagator will (at most) propagate a single 
two particle bound state and a composite two particle state. The most general ansatz with momentum conservation 
consistent with this physical picture is given by: 

A{j>x,P2,P3,Pi) = $(pi -P4)6(p 2 -ps)F{p l ,p 2 ) +6(p! +p 3 -p 2 -p4)G(p 1 ,p 2 ,p 3 ,p 4 ) 



r< + * 



(39) 

Inserting this into (|3"§|), we immediatly find: 

F( Pl ,p 2 ) = a ( Pl )a (p 2 ) 

g 2 1 g 2 f d d k 

G(pi,p 2 ,P3,Pi) = -i^j^y^(PiW°(P2)cr (p 3 )a (p4) - i—a°(p 1 )a°(p 2 ) J j^-^G{k,p 2 - pi + k,p 3 ,p 4 ) (40) 



Iterating this second equation for G(pi,p 2 ,p 3 ,pi), it is not hard to see that we are reproducing the series expansion 



for: 



2 

G(pi,P2,P3,Pi) = -«^-cr (pi)o- (p 2 )cr (p 3 )cr (p 4 ) 



1 



l + iH^ (k)^(k+p 2 - Pl ) 



(41) 
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It is easy to understand this last equation in terms of more familiar approaches to the large N limit: The object 
J a a is a bosonic bubble. Thus G(pi,p2,P3,Pi) is simply a sum over chains of bubble diagrams. The factor in square 
braces is in fact the propagator of the auxiliary field usually introduced to study this model 0. Notice however, that 
our a propagator A(p\,p2,P3,Pi) consists of two terms. The first term, which has no analog in the auxiliary field 
approach, is crucial to obtain a systematic expansion. 

The mass spectrum of the theory is determined by searching for the poles of the propagator A(pi,p2,P3,pi). 
Rather than performing a full non-perturbative analysis of the spectrum, we content ourselves with constructing 
A(pi,P2,P3,Pi) to 0(g 4 ), by iteration. A straightforward calculation shows that the connected piece of A(pi,p2,pz,Pi), 
to 0(g 4 ) is given by 

Aconnf \ 9 2 S (Pl ~P2 + P3 ~ Pi) Qf \ 0/ \ 0/ \ / \ , ( -\2 9^ $(Pl ~ P2 + P3 ~ Pi) 

A conn { Pl ,p2,p3,Pi)= -i-, (T u (pi)a u (p2)cr u (p 3 )<J (Pi) + (l) 



(2ir) d " yl ' 1J " yl " ±/ ' vv 16 (2tt) 

a (p 1 )a°(p2)a (p 3 )a (p i ) J + p 2 Px)<r°(k). (42) 

Inserting the expression for the leading order (|34|), yields complete agreement with the 0(N) term in equation (|3C|), 
as it should. 



5. Corrections to the one and two point functions of the effective theory. 

In order to compute the 0(— =) correction to the effective theory one point function, we need to compute the 
tadpole, and cubic interaction vertices of the effective field theory. These interactions arise from the actions 

St, = iLV{0)^=B s 



8. = -« ' * (43) 



respectivley. For completeness, we also present the j? quartic and subleading quadratic vertices 

S sq = -iL d 6 d (0)^-B 2 



1 

'4/V 

Standard techniques yield the following Feynman rules: 



= *—Bi (44) 



L d S d (0) r j-li 
Pl,P2 VN [ 



kl,k 2 Pl,P2 LV(0) r o-ln r^O-li , 



^(K _1 ] fe2Pl W°-\ 2qi [^\ 2kl + [a -\ 2kl [cr*-% 2qi [<J G ~\ 2P1 ) 



3lVW 



ki,k 2 Pi,P2 
9i,92 h,k 4Dv 



-MK^W, [^-\ 2qi [a°-\ 2h [a°-\ kl + [a -% 2Pl [cr " 1 ]^ [a ' 1 }^ [a " 1 ]^ 



ki,k 2 Pi,P2 



+ ]k 2qi [v° \ 2 iA^ }hpt[<r° + X^A^ l ] q2Pl [(r° ] P2 hW° 1 ]; 2 fe 1 

W°~\ qi [* 0_1 W W°~X 2kl + [o ' 1 ]^ W°^}i 2Pl W°-\ 2qi [a Q -\ 2kl ) 



(45) 



In the above, we have employed an obvious matrix notation. Recall that a°(k,p) is diagonal in momentum space, 
so that 

o-i _ 5{k-p) 
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Now, we turn to the computation of the (3(-A=) correction to the one point function < a >. Two processes 
contribute: the cubic tadpole and the linear tadpole. Using the expression ([|2]) for the propagator, and the above 
rules for the vertices, we obtain the following expressions for the cubic tadpole 

o 

-^fdpa°-\p)A(p,p, PuP2 ) + 5(p 1 -p 2 ) J ^a°(p 1 )a°(p 2 )[f 1 ^a°(k) 

d d k „OfU\ r d d l „0fl i U m \„0/J\ ig 2 f d d k „Q ( h\ f d d l ( „0 



(47) 



up to 0(g A ) and for the linear tadpole 



1 



^-jMjdpa" 1 (p)A(p,p,pi,p 2 ) 



(48) 



Notice that separately both the linear tadpole and cubic tadpole contain momentum dependent infinities. The sum 
however is a well defined quantity, with all divergences proportional to L d S d (0) cancelling and 



■ 2 

< 77 >= 5(pi -p 2 )^|=^ (piV (p2) 



(2tt)<* r {k) 4 



d d l 
(2V^ 



{a°(l)) 2 a°(k) + a°(k)a°{l)a°(k + l- Pl ) 



(49) 



Inserting the expression ( |34| ) for cr°, and recalling that our fields are rescaled by a factor of N, one obtains complete 
agreement with the perturbative result (p8|). This analysis shows that the linear tadpole is essential to obtain 
agreement with the Feynman perturbative analysis and to remove momentum dependent infinities. 



III. FERMIONIC VECTOR MODELS. 



In this section we show that the approach followed in section 2 for bosonic vector models also applies to fermionic 
vector models. This is an important observation and this is not only because fermionic systems have more important 
physical applications as we now explain. 

By and large there have been two approaches to the large N limit of vector models: in the first one, auxiliary 
fields are used. We have already mentioned in the introduction some of the shortcomings of this approach. The other 
approach is based on the collective field theory Q| idea of changing variables to invariant quantities. However as it has 
been discussed in section 2b the equation satisfied by the appropriate Jacobian ultimately stems from a Hermiticity 
requirement of the Hamiltonian. It is not obvious how to generalize this Hamiltonian approach to fermionic systems, 
although some partial success has been achieved in terms of pseudospin variables ||] . We will not pursue this method 
in this article but will show that, provided that one is prepared to consider time ordered product expectation values, 
the fact that the Schwinger-Dyson (loop) equations of the theory imply a differential equation for the Jacobian, as it 
was demonstrated in section 2a, is straightforwardly generalizable to fermionic systems. Moreover our approach will 
be fully covariant. 



A. Schwinger Dyson equations and Jacobian 



We assume again that we are dealing with U(N) invariant actions i.e. actions invariant under ip® — > ip' a a — U ah i\) a , 
-> -0a° = U ab ip^ b with U ab an arbitrary U(N) element, a is a Dirac index. The invariants are now 

a aP {x,y)=r a {xW {y) (50) 
and as before we are interested in time ordered product expectation values of invariant opertators 

< F[a] >=< [[a ai0i ( Xi , yi ) >= Jd^d^s ( 51 ) 
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A set of Schwinger-Dyson equations for the invariants follow from the identity 

6 



= / Di/j*Dip 



5^{x) 



r p (y)F[a]e^ (52) 



It is important to remember that the fields ip have to be treated as Grassman variables in all manipulations to 
follow. We obtain 

< N6 Pp S d (x - y)F[a] > - < > - < r p (v)^F[a] >= (53) 

Postulating the existence of a Jacobian J defined by 

; ■ • / ' / a '/ s = I Da. J F[a]e lS (54) 
we consider the identity 

0= f Da f d d z—^- -(a ap (z,y)JF[a]e iS ) (55) 

J J oa a p(z,x) 



which implies 



< m6 0p 6 d (x - y)5 d {Q)L d F[a] + I d d za ap (z,y) 9lnJ F[a] + f d d za ap (z,y)- '' FM 



da a p(z,x) J ' da al3 (z,x) 



+i I d d za ap (z,y) a _ F[a] >= 0. (56) 



dS 

' da a p{z,x) 

In the above equation m is the dimension of the representation of the Clifford algebra. Using the chain rule 



5ip%{x) J J S^(x) Sap p (y,z) J p 8a 0a (y,x) 
in equation (153) it becomes 



/& F f S S 
dza ap (z,y)- > -i < dza ap (z,y)- -F[a] >= 0. (58) 

ba a p{z,x) J da a {3{z, X ) 

Comparing this last equation with equation (^6|) for arbitrary F[a] we obtain 

Sin T 

d d za ap (z, y)- = -(N + mL d 6 d (0))6 Pp S d (x - y), (59) 

6a a p(z,x) 

The solution to this equation is 

InJ = — (JV + mL d 5 d {Q))trTrlna (60) 

The trace in the above equation runs over both Dirac and functional spaces. This is the main result of this article. 
As mentioned earlier the leading term of this Jacobian has been obtained in 0. This Jacobian should be compared 
to the bosonic one Jl4|). 

It is known that in models such as the Gross Neveu model an enlarged 0(2 N) symmetry is present that is better 
exhibited in terms of Majorana components |Tsfl . One could have considered a set of invariants 

a' ap {x,y)=r a {x)my) (61) 

for which the Jacobian is easily seen to be 

InJ = —(N + mL d S d (0))trTrlna (62) 

Finaly if flavor degrees of freedom are added as in Nambu Jona-Lasinio type models then one needs only include a 
further trace in flavor space in the Jacobian (pQ). 
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B. Effective action and large N configuration. 



Our effective action is given by 

S eff = -ilnJ + NS = iNtrTrlna + NS + imL d 6 d (0)Trlna = NS + Si. (63) 

We have assumed that we have again rescaled the fields so that an overall factor of N multiplies the action. The 
leading N — > oo configuration is determined by 

fW0 (64) 

Shifting about a as: 

cr a p(x, y) = ala(x, y) + -^=T] al3 (x, y) (65) 
v iv 



we expand 



S eff = NS (a°) + - l -C 2 + \D 2 + | [^l)" +1 [^y + ^(0)^1 + (66) 



where 



and 



D„ = / d d xi... I d d x n I d d yi... I d d y n - ^-f -, s\u Q iloL- L 1 {xi,yi)...'n oln p n {x n ,y n ) (67) 

oa ai f3 1 {x 1 ,y 1 )...da an p n (x n ,y n ) 



C n = trTr^Y^T (68) 



At this point, a few comments are in order. The original fields ip and ip are fermionic fields, and consequently they 
become Grassman variables under path integral quantization. Using the Grassman nature of the original variables, it 
follows that 7ct 0! (T o-/3 is an antihcrmitian bosonic bilocal; all knowledge of the original fermionic statistics is now coded 
in this property and the specific form of the interactions in (|66|). This provides a non trivial check of the effective 
field theory that can be carried out with perturbation theory. 

1. Perturbative Results 
We specialize now to the following Lagrangian density (written before rescaling): 

£=*(^-m)$ + ^(t*) J + |(M a (69) 

For now A and A5, assumed to be of the same order, are left arbitrary, although later the cases A5 = (the Gross 
Neveu model) and the case A 2 = —A 2 (a Nambu-Jona-Lasinio type model) will be considered in detail. 
By standard diagrammatic techniques, one can obtain for the two point function 

< m*)Mv) >= J ■^3 G ^( k >~ ik{x ~ v) ( 7 °) 

with 

G^ a0 { P ) = N( + Q + Q Q + 8 +0(A 6 )) 

+( O. + OQ + OO + . 



+ 



-0(A 6 )) + 0(i 7 ) 



(71) 
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We have introduced a dot into our notation, so as to indicate the contraction of Dirac and color indices. Since the 
summation over color indices is taken into account in the overall factor of N, each dot should be thought of as a sum 
over the Dirac matrices 1 and 7 s , multiplied by the interaction strengths A and A5 respectivly. For instance 

Q = Sm=A1,A 5 7 5 T^m. I jh?)Z Tr \- M T^l\ M Y^rn 

(72) 



For completeness, every single diagram appearing in ( |7l| ) is explicitly written down in appendix A. 
The connected piece of the four point function is given by: 

afc.-™,*,.™,,)^)™^^^^, (73) 

where 

G^ COnn {pi,p 2 ,P3,P4,) = N(/\ + (•) + X. + permutations + 0(A 6 )) + 0(1). 



Explicitly: 



-Q— — (*) 2 Z)m,JV=A1,A 5 7 5 [ fi-m^ tf 3 -m\p' a ' I (Sir) 1 * "^M f6-m^ fl-fa+fa-m^^ fa-m^ fa-m\p a 



2. The Leading Configuration 

With a translationally invariant ansatz: 



we obtain from 



where 



a = 



(2n) d ' J (2tt) 



(74) 



(75) 



(76) 



(77) 



a " '»''''■ n ' = I -0^e ik ^a a0 (k) (78) 



iv° P a(k) = (fi-m + X 2 a + A^cfs)^ (79) 



'^iaUk) = -G = - I ^L T r[G^{k)] (80) 



and 
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05 



d d k 



(2tt) 



d'y a /3 cr af3(k) = -G 5 



d d k 



(2tt) 



^Tr[ 7 5 G«(fc)] 



(81) 



This is the standard gap equation. We again assume that all integrals are regularized. Using (|79|), we can write, to 
0(A 4 ) 



i-^(\ 2 G + A 2 7 5 G 5 )-^- + W 2 -^(A 2 G + A 2 7 5 G 5 )-^(A 2 G + A 2 7 5 G 5 )-^ 

m p — m p — m p — m p — m p — m 



(82) 



Since 



A 2 G + A 2 7 5 G 5 = 



dfik i n v-^ f d d k f d d k' 



(2?r) d L #-m J 



M,AT=A1,A 5 7 5 



5 7 (27r) d 7 (2?r) d 



-]Tr[JV- ]M + 0(A 6 ), 



m ft — m ft' — m 

it is straightforward to substitute this expression into (182) to obtain 



(83) 



nil 



-0(A 6 



(84) 

in agreement with the leading order term in the diagrammatic expansion for the two point function. Notice that in 
this case, not all diagrams contribute with the same sign, due to the appearance of closed fermion loops in the leading 
order. It is remarkable that the original fermion statistics is completely accounted for by the change in sign of the 
leading order term in the Jacobian by comparison to the bosonic case. Recall that 7° r CT Q , ( g is an antihcrmitian bilocal 
bosonic field. We will see that this change of sign correctly reproduces the Fermi statistics to all orders of the l/N 
expansion. 



3. Effective Field Theory Spectrum. 

In this section, we use the leading order ff^, at the level of quadratic fluctuations, to obtain the mass spectrum of 
the effective field theory || . From (|6rj ) , the leading quadratic action is 

S 2 = - l -C 2 + ^D 2 (85) 
= ^ J dxidx 2 dx 3 dx i [(\ 2 S pa 6 p > a i + \lj pa l p i a >)o~(xi - x 2 )6(x 1 - x 3 )S(xi - x A ) 



■la 



-" (xi,x 2 )] ap >[a° ± (x 3l x i )] a ^ p ]ri pa (x4 : ,xi)ri^ a ,(x 2 ,x 3 ) 



Now, using a symmetric Fourier transform 

Vp«(*i,"») = J ^^^^Vbi^) (86) 

the quadratic action is 

1 f A 2 A 2 

S 2 = ^ J dpidp 2 dp 3 dp i [j^^5 pa 8 p , a ,8{pi+p 3 -p 2 -p i ) + j^yl b pa ll' a ' 5 {Pi + P3 - Vi - Pi) (87) 

-iS(pi -p 2 )[<T° (pi)]ap'5(p 3 -P4)W° 1 (P3)} a 'p]Vpa(P4,Pl)Vp'a'(p 2 ,P3) 

in momentum space. Using this action we find that the propagator A apa i p i of the field a ap satisfies: 
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'd Pl d P2 [(\% v 6 pa + *hlull / {kl+P ^y P2) ~iS(h P2 W va 1 (p 1 )6(p 1 - k 2 )a° pfl \p 2 ) 



Using the argument preceding (39), we make the ansatz 

A p y pT (pi,p 2 ,P3,P4) = 5{pi -Pi)8(p 2 ~ P3)F f _ l v pT (pi,P2) + 5{pi + p 3 ~ Pi - Pi)G p V pT {p\ , P2,P3,Pi) 

= rx + * 

Inserting this ansatz into (|8^), we find the solution 

(27T)° 



(89) 



A a >p>ap(p3,P4,,Pl,P2) = S(P3 -P2)S(j>l ~ Pi)°a> p{P2)o° apl (pi) + JTT^d^ipi ~ & + P3 ~ Pi) x 



where 



and 



X J2 {Pi)Ma° 1 {p3)]p^'T MN {p 2 -p 1 )[a 0l {p 2 )Na 0l {p 1 )]p a (90) 

M,JV=A1,A E 7 5 



F M n(p2 ~Pi) = [1 + iG(p 2 -Pi)]mn ( 91 ) 



Gmn(P2-Pi) = J -^Tr[M<j oT (k+p 2 -p 1 )Na QT (k)]=G NM ( P i~P2) (92) 

This last term corresponds to summing up the bubble diagrams, which we know dominate the large N limit. The 
diagonal elements of Tm n are identical to the propagators of the auxiliary fields which are usually employed to study 
this model |Q. The first term, consists of disconnected diagrams in terms of the original fermion fields and has no 
natural counterpart in the auxiliary field approach. This term is however, crucial to obtain a systematic expansion. 
This will already be evident at the perturbative level, in the discussion of the next section. It is easy to see how the 
Fermi statistics is reflected in the above propagator: The first term corresponds to an exchange of the two ingoing 
fermions, with respect to the leading configuration. It thus appears with the opposite sign to the corresponding term 
in the bosonic propagator (|39|). Also, the fermionic bubble in the denominator of the second term appears with a 
minus sign, reflecting the fact that closed fermion loops come with a factor of —1. 

A full nonperturbative treatment of the spectrum will be discussed in a later section. At this point we simply 
construct the connected piece of A a i p i ap to 0(\ i ). First, expand the second (connected) term of ( po| ) as 

ATp^pi p 3 ,Pi,Pi,P2) = jAdS(pi-P2+P3-Pi) E lv oT (PiWa oT (p3)}p, a '{v oT (P2)M* oT ( Pl )}p a (93) 

A/=A1,A 5 7 B 

- MjSiPi -P2+P3~ Pa) Yl \v QT (Pi)Ma QT {P3)] P - a -G M N(P2 - Pl )[a oT (p 2 )Na oT ( Pl )}p a + 0(A 6 ). 

^ ' A/,A r =Al,A 5 7 5 

Using the expansion for (<x°)J Q = er^ = — G^l to 0(A 2 ) obtained in a previous section, one readily verifies that 
the above expression is in exact agreement with the expansion j7^-[T7|). 



4- Corrections to the one and two point functions. 

In this section we compute the tadpole, (subleading) quadratic, cubic and quartic interaction vertices of the effective 
field theories. We use these vertices in this section to obtain agreement with the perturbative results quoted in 
subsection 2 for the (fermion) two point function. These vertices will be used in a calculation in a later section. 

From the following terms in the action 
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S tp = imL d 5 d (0)-^Ci 
V N 



1 



3ViV 



C 3 



S s , = -*mZ d S d (Q) — C 2 



2N 



S a 



(94) 



we easily derive the rules 



Ba,0 



Pt,P2 

a/3 jis 
k\,k 2 ^ pi,p 2 



mL d S d (Q) r O" 1 ! 

^rjj V 7 \p2Pp1a 

m L d 5 d (0) r Q-l 



2N 



f O- 1 ! 



a/3 



a/3 7^ 

21,92 ^1,^2 



</>-0 / \//T 
^1,^2 Pl,P2 



3!\/]V^ Cr0 Ife^Pi^t' 70 ]p2-09l Q [ <7 ° ]<32/3fei7 + I " ]p2"0' £ l7 1 170 ]fc2^<Jia I* 7 * 3 ]<?2/3pi^) 

4T/v([ <7 ° ]fc2i/ , PlA 1 [ Cr ° ]p2T<jia[ ,T0 ]92/3/i7[ <7 ° ]/2ffel0+[ cr ° ] fc 2 0Pl M I ' ]p2"r/i7 X 
X[(T° ]i 2 i/g lQ [o-° 1 ]q 2 l3k 1 <f> + [<7° ]k 2 ^ qi a[o° ]q a /3h-/[<7° ]hvpm[<T° Wfci0+ 
+ [ cr ° 1 ]k 2 iPq 1 a[v° 1 ]q 2 /3pi^[ cr0 1 ]p 2 Tl 1 j[v° \hvknf, + 1 ]/c 2 ^/i7[ (j0 ta^gia* 

X [ cr ° ]g2/3piM[ CT ° ]p2rfci0 + [f° Ifea-i/'ilTl^ ]hvpitt\- a ° }p2Tq 1 a[ cr0 ]q- 2 f3k l ^) 



We have again gone over to an obvious matrix notation. Since a is diagonal in momentum space, we have 

[c° ]pakj3 = S(p - k)[(J° 1 (p)] a /3 



(95) 



(96) 



As in the bosonic case, two processes contribute to the 0(-^=) correction to the one point function: the cubic 
tadpole and the linear tadpole. Using the expression (|9Cj) for the propagator, and using the expansion (93) (for the 
second connected term of the propagator), together with the above rules, we write down expressions for the cubic 
tadpole 



o 



= mL % ia) S d d pa Q p e 1 {p)A epafi {p,p, hi, k 2 ) + -^5{k x - k 2 ) [t / jgfa Ew=ai,a»v 
[a° T (k 2 )Ma oT (p)Ma° T (k 1 )} l3a + J2m,n=Xi,X 5 ^ J j^[a° T (k 2 )Na° T (p)Ma° T (h)} 0a x 
G M N{k 2 -p)+fj$frj ^Tr[Ma oT {p 1 )Na° T (p 1 )Ma oT {p 2 )}[a oT (k 2 )Na oT (k 1 )} 0a 



(97) 



to 0(A 6 ) and for the linear tadpole 



1 



mL d S d (0) p Jd__0-1 



J d d pcr Q pe (p)A epa/3 (p, p, h,k 2 ) 



(98) 



Again we find that the linear tadpole is essential to cancel momentum dependent infinities arising from the inte- 
gration in the cubic loop. The sum of the two tadpoles is given by 
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(2tt) 



< Va p(h,k 2 )>=4=S(k 1 -k 2 )[i [ £ [a oT (k 2 )Ma oT (p)Ma oT (k 1 )} f3a + £ f 

x[a oT (k 2 )Na oT (p)Ma° T (k 1 )} (3a G MN (k 2 -p) + J J ^Tr[Na° T ^Ma * (p^Na * (p 2 )] x 

x [a 01 \k 2 )Ma° 7 \h)} 0a + 0(X 6 ) (99) 

Substituting the expressions for tx° and Gmn given in (|8^ ) and (|9^), we find complete agreement with the O(l) term 
of (0). Notice that once again all diagrams enter with the correct signs - thus corrections computed using our Jacobian 
( |60| ) have Fermi statistics correctly accounted for! (It should be recalled that < a a p(k\, k 2 ) >— —G^?(ki,k 2 )). 

IV. A NON PERTURB ATIVE ANALYSIS OF THE EFFECTIVE FIELD THEORY SPECTRUM 

In this section, we obtain the non perturbative effective field theory spectrum and propagator for the Gross- 
Neveu model and a Nambu-Jona-Lasinio type model (a model with a continuos chiral symmetry) and relate them to 
previously obtained results. 

A. The Gross-Neveu Model 

For the Gross-Neveu model, we set A5 = m = and work in d = 1 + 1 dimensions 

Notice that this model has a discrete chiral symmetry 

< -» ilM r a ^-Vuil a (101) 

Using the formalism developed in section we find that the leading configuration a satisfies 

-i[a°~X a {k) = {^ + X 2 a) ap (102) 

where 

^Tr[a{k)] (103) 
Inserting ( |102[ ) into ( |l03| ) we find the standard result 

where A is an ultraviolet cut off. Therefore 

\ 2 a = Ae-f? (105) 



Thus the fermions have aquired a dynamically generated mass X 2 a and the discrete chiral symmetry ( |l0l| ) has been 
spontaneously broken, as is well known Q. This mass can be made finite as A — ► 00 by adding the counter term 

Set = ~ [ -^^r^-d 2 p 3 d 2 p 4: S(pi + p 3 -p 2 - P4,)rj aa (pi,p 2 )r}w(p3,P4:) (106) 



2 J 2?r 2tt 

to our effective action, corresponding to a coupling constant renormalization of the original theory. Fixing the 
renormalized mass to some value m r 

m r = Ae~^?~c = Ae (107) 

we can fix c 
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1 m 2 

c = A * + [±log!£.]- 1 (108) 
Using the arguments constructed in section B3 we know that the effective field theory propagator 

\ V ^'u'{Pl,P2,P3,P4,) =< V^{Pl,P2)V^'u'(P3,P4) > (109) 

can be expressed in terms of the leading configuration <t° as 

\ V ^ V '(P1,P2 > P3P4,) = S(p 3 ~ p 2 )S{p 1 -p4)o-^(pi)<T°/ J/ (p 3 ) 



— ^$(pi+P3-p2-pi) — w jr M — (HO) 

(2?r)2 ^ P±> l+iX 2J^^ {k)(J Ak _ pi+n) ^ > 



To construct the propagator explicitly, we perform the integral 

d 2 l 
(2n) 2 



A 2 , /4m2_fc2 TwTp _ ^Zfc2 An, A , , 



where k = pi — p 2 - Now, from (|107[) we see that 



so that 



^log — = 1 (112) 

7T 771 r 



2ttV -fc 2 ° ,9 V4m 2 -/c 2 + V^F J 



Inserting this into the propagator ( 110 ), we obtain 

1 2 

A^'u'ipi)P2,P3Pi) = S(P2 ~ P3)5(P4 ~ Pl)[ 



<5(Pl +P3 ~P2 -Pi) 



2ir 1 $2 — Tn r p'l — m r v> * p 1 ^ — m r ^3 — m r v M J 4mj-t ! , ry/4ml-k 2 -V-k 2 - 



(H3) 



(114) 



Notice that all dependence on A has disappeared as it must. 

The interpretation of the above propagator is clearer in terms of rapidity variables p° — mcosh(6)iP l — msinh(6) . 
If <f) = Q\ — O2 the connected piece of the propagator (after the removal of the external legs) is given by 

U(0) = tJ^P-. (115) 

For on-shell particle-antiparticle scattering (corresponding to time flowing from right to left in (89)) we let <fi — * <fi—iir 
(pi — * pi,P2 — ► — P2) and obtain the amplitude 

, N coth($r) 

DU) = 2m- — (116) 

(j) — 17T 

Since s — (pi + P2) 2 = 4:in 2 cosh 2 (^) we see that as 4> — ► the above amplitude displays a cut at s = Am 2 . 



17 



B. A Nambu-Jona-Lasinio type model 



The model we consider in this section has A§ = — A 2 and m = 

l = ri @r + ^{rr? - ^(ri 5 r) 2 (m 

We again work in d — 1 + 1 dimensions. This model is invariant under the global continuos chiral transformation 

C^e^C r a ^-i> a v e a ~<\ a (118) 

The leading configuration a satisfies 

- i[<r°^U(k) = (# + A 2 a + A 2 7 5 a) Qp (119) 

where 

<hk ;Tr[a{k)] (120) 



(2ny 
and 

d 2 k 



(2^) 2 

Inserting ( |l 1 9| ) into ( |l20|) yields the standard gap equation p9| 



Tr[7 5 a(fc)] (121) 



[^J + AV^ + AV]- 1 = / -_[_____-] (122) 



(2tt) 2L ' 1,1 ivu J (27r) 2L fc 2 -4A 4 (a 2 + a 2 ) 
Similarily, the gap equation for a is obtained by inserting ( |ll9| ) into (|l2l| ) 

rf 2 fc r -»7 5 1 / ^ , 4»A 2 a 

(27r) 2[ /(; + A 2 7 5 a + A 2 CT J,y,y J (2tt) 2 [ £; 2 - 4A 4 (<7 2 + ct 2 ) J 1 J 

Fixing the mass to some value m r , and parametrizing [fl9| 

2A 2 er = m r cos9 2X 2 a = m r sin6, (124) 

we find 

1 1 A 2 , , 

2A^ = 2^ (125) 

The expression for the effective field theory propagator in terms of the leading configuration is 

+ j3- S{pi+k - P2 - P4) v ^fa^fwwi^y^faw (126) 

(2,P „,£i, 7 . 1 + i ^ /^rr[^ T (*)M^ T (*-j> 1+ »)W] 

To make the following arguments simple and transparent, we pick 6 = in cr^g- To obtain an explicit expression 
for the propagator, we need to compute three different integrals. The first integral was considered in the previous 
section, so we simply quote the result Q 

/d 2 k 
—-~2(T^(k)a y(T (k - pi +p 2 ) 



2tt V -k 1 L ^4m 2 - fc 2 + V-fc 2 
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The second integral 



h = 



<hk Tr[a oT (k-p 1 +p 2 )la° T (k) 1 5 } 



(27T) 



(128) 



vanishes, due to the trace over Dirac indices. The third integral is easily computed 

d 2 k 



(27T) 



r Tr[a° T (k - p x + p 2 )ij 5 a° T (k)ij 5 } 



A 2 Um 2 -k 2 1 r V4m2-P- V^F n Ag . A 2 A 2 



' ^J4m 2 - k 2 + V-k 2 27r m 2 



(129) 



where A is an ultraviolet cut off. Using the renormalization condition, in the form (125), we may rewrite this as 

(130) 



An 4ml 



-log 



\J Aui 2 - k 2 - -J-k 2 



2ttV -k 2 L y/4m 2 -k 2 + V^ 
Thus, the effective field theory propagator reads explicitly 

A i1u ^ v >(pi,P2,P?,Pa) = S(p2 ~P3)S(pa -p%)[- 



ii — m r p3 — m r 



S(pi +P3-P2 - Pa) 



2TT $2 — m r i>\ — m r p'i ~ m r i>Z ~ m r " ** , / 4mjS-fc 3 , f ^mg-t 2 - 



-k 2 - 



-k 2 



8(Pl +P3 ~P2 - Pi) 



2n $2 — r^r 



I fir V> * $4 — m r i>z — m r " M / 4m 2 -k 2 ^^r ^/4m^-fc 2 -\/-fc 2 



(131) 



-k 2 



V 



4 m 7 



-k 2 ' 



The first term again represents the crossed propagation of two free fermions. The second term has been discussed 
earlier. The third term above, is however new. It has the property that 



4m 2 . — k 2 
~T 2 



log[ 



V^k 2 - 



yjAm 2 - k 2 + V-fc 2 ' 



=n=2. 



\k 2 =a 



(132) 



This has been interpreted as the signature of a massless scalar in the spectrum of the model pQ ]. This massless 
particle is associated with the fact that the continuos chiral symmetry of the model, is "nearly broken" |21|. 
Finally, we remark that the Lagrangian (117) with added flavor degrees of freedom 



c = jhl> + -a 2 [(^/0 2 - $7 6 r,vO(fy , W)] 



(133) 



is a popular candidate to study the low energy phenomenology of the light mesons. In the above, the Tj are taken 
as the generators of SU (2) for the two flavor model, and as the generators of SU (3) for the three flavor model. The 
invariant correlators now carry four labels 



*%{x,y)=rj{x)i,l i {y) 
where i,j are flavor indices. It is not hard to construct the leading configuration 



d 2 k 



-ik(x-y) I 



1 a/3 



where the mass m 2 = A (a + a ) satisfies the gap equation 



1 



4zA 2 



d 2 k 



(2tt) 2 k 2 - 2\ 2 m 2 

The effective field theory propagator is easily written in terms of this leading configuration 



(134) 



(135) 



(136) 
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tAg ^ [T 0r (Pl)A^ 0r (P2)]^[T 0T (P3)^ 0r (P4)]g, 

(2?r)2 + p 4 ) 2. ! + ^ J ^ Trtr[a 0T Ma0 T {k _ n + p2)N] 



where Tr denotes a trace in Dirac space, tr denotes a trace in flavor space and M, N are now direct products of 
matrices belonging to Dirac space with matrices belonging to flavor space. fj] Notice that in this case a triplet of 
massless pseudoscalar bosons appears in the spectrum. These are usually interpreted as the 7r + , the ir~ and the ir° 
particles. 

In the context of this phenomenological model, our effective field theory is nothing but the theory of the mesons 
built from quark anti-quark pairs. The couplings of the various meson-meson interactions are proportional to N (the 
number of quarks) raised to some negative (integer or half integer) power. 

C. Homogeneous Bethe Salpeter Equation. 

We have shown how to compute the Feynman rules and propagator for our effective field theory and have verified 
their correctness perturbatively, and non perturbativley to first nontrivial order in . However before any scattering 
amplitudes can be computed, we have to supply a set of asymptotic states. S matrix elements are then taken using 
these asymptotic states, as usual. 

In any field theory, only the quadratic term in the action provides a harmonic oscillator and thus a spectrum 
consistent with a set of free particles. Thus, it is the quadratic term that codes the asymptotic states of the theory, as 
solutions of the corresponding homogeneous equation. This should always correspond to the large N approximation 
to the homogenous Bethe Salpeter equation. 

For example, the effective field theory wave functions for the Gross Neveu model satisfy 

A 2 f 

i(i>2 - m) p > a ri pa (pz,p2){i>z - m) pa < + J dp 1 dp 4 S(p 1 + p 3 - p 2 - Pi)il aa (Pi,Pi)5 a ' P ' = 0. (138) 
It is possible to obtain a particular solution to this equation. Making the ansatz 

Vpa(P3,P2) = Vpa(P3,P3-P2), (139) 

we find that fj satisfies the equation 

A 2 f 

(02 - m) p , a fj pa (p 2 ,P3 ~ P2)(p'3 ~ m)pa> =i ^2~ / d PlVaa {Pi , P3 ~ P2)5 a > p> ■ (140) 

Rewriting this last equation as 

i i A 2 f 

fj^(p2,P3-P2)=[-, r-: L„*TT / d PiV aa (Pi,P3 ~ Pa) (141) 

p2 — m) pz — m J f v vk 1 J 

taking the trace of both sides and integrating over p 2 keeping p% — p2 = k fixed, we are lead to the consistency 
condition 

1 = *-T% / dp 2 \-r^ 3 1 1 (142) 

This condition requires k 2 = 4m 2 which, since k is the physical momentum transfer for the particle antiparticle 
channel, corresponds to the leading order mass shell condition for the fermion-antifermion bound state. 



There is a slight abuse of notation in (137) - please refer to (90-92) for clarification 



20 



V. SUBLEADING i CORRECTIONS. 

The formalism which was developed in sections II and III can be used to obtain sublcading corrections to any 
correlator of interest. For instance the correction to the propagator results from the following diagrams in the 
effective field theory 



< SaSa >= 



1 



(143) 



with the Feynman rules given by (95) in terms of the leading configuration <to and where the propagators in the 
diagrams refer to (114). It is essential that both the disconnected and connected pieces of this propagator are 
included. One observes that the contributions from the subleading term of the Jacobian (proportional to L d S d (0)) 
precisely cancel similar infinite contributions generated in the other diagrams. All the remaining expressions then have 
an interpretation in terms of original Feynman diagrams. In this precise sense the subleading term in the Jacobian 
provides a normal ordering. Diagrammatically, we find: 






+ 




(144) 

In the above, we have only included basic skeletons; permutations of these basic skeletons need to be included. For 
example, the first diagram appearing would be repeated with another diagram in which the second fermion is dressed. 

It is well known that the Gross Neveu model is exactly integrable with no particle production (although it has a 
rich spectrum ]l8|]) and the exact S matrices are known mM . For two particle scattering the S matrix element is given 
by 



< Pb(Pl)Pd(p2)\Pa(pl)Pc(p2) > l 



S bd {e)5{f,\ - P \)5(p\ - p\) - ac S db (6)5(p{ - p\)5(pl - P \) 



(145) 



where 
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acS bd (6) = a 2 (6, 2N)6 ab S cd + a 3 (d, 2N)5 ad 5 bc , 



(146) 



In the above equations 9 = &\ — #2 where Q\ and 62 are the rapidity variables of the incoming particles. The 
above U(N) symmetric two particle S matrix is consistent with the underlying 0(2N) symmetry of the Gross Neveu 
model (22]. 02 and a 3 are the standard symbols used to describe the exact S matrix with O(N) symmetry. 

It turns out that the simplest way to verify the correction to the propagator is to consider the contributions 



of the diagrams in (144) to forward two particle scattering corresponding to time flowing from top to bottom. The 

ab S ab (forward) = N 2 (a 2 + ^-a 3 ). (147) 



corresponding S matrix is then 

1 



N 



By looking at the diagrams in (144) we see that the first two diagrams can not contribute to forward scattering. 
The third diagram is clearly the leading (O(jj)) contribution to 0-3. The remaining diagrams should sum up to the 
contribution to 02. They are precisely the diagrams considered by Berg et. al. ]l3| who have indeed confirmed 
that this is the case. 



VI. CONCLUSIONS 



In this paper, we have used colorless bilocal fields to study the large N limit of both fermionic and bosonic vector 
models. By requiring that the Schwinger Dyson equations in terms of the original variables agree with the equations 
derived directly in terms of the invariant variables led to a functional differential equation for the Jacobian. The 
equation was solved exactly, leading to an exact effective action. This effective action was then shown to reproduce 
the familiar perturbative expansion for the two and four point functions. In particular, in the case of fermionic vector 
models, the effective action correctly accounts for the Fermi statistics. The theory was then studied non-perturbatively. 
The stationary points of the effective action provide the usual large N gap equations. The homogeneous equation 
associated with the quadratic (in the bilocals) action is simply a two particle Bethe Salpeter equation. Finally, the 
leading correction in i was shown to be in agreement with the exact S matrix of the model. 

There are a number of interesting questions which can now be asked. Firstly, it is clear that the above invariant 
variables are classicaly commuting functions. This change of variables therefore provides a bosonization valid in an 
arbitrary number of dimensions. It is interesting to ask if a link can be drawn to more conventional bosonization 
schemes. Also, this bosonization may prove to be a powerful tool for analyzing many body condensed matter systems. 
Secondly, it is well known that the large N limit of the Gross-Neveu model posseses a number of interesting config- 
urations, corresponding to many particle bound states. The present formalism is particularly well suited to studying 
two particle bound states. It would thus be interesting to determine the Feynman rules of the effective theory in a 
two particle bound state background. One would expect that these rules may be simpler than the rules obtained in 
this paper. Finally, the possibility of obtaining an equal time approach for fermionic vector models remains of great 
interest. 
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Appendix A: Diagrams to integrals. 

In this appendix, explicit expressions are supplied for the fermionic diagrams appearing in the text (see equa- 
tion@): 



(148) 



=-iEM=x,^^I^Tr[M J ^]M T : 



(149) 



(150) 



(151) 



O — iV i r d d k \/r i u i_ 



(152) 



(153) 



(154) 



~7 2 V i f d d fc 1| i AT i r d d k' J, rjy T j ]» t 

' Z^M,iY=l,<7 5 j/- m J (27i-) dJVJ ^-m JV jC-m J (2^)^ ' ^ JV JC-mJ^ j/-m 



(155) 



(156) 



;2V S il/f f ti d fc d d k' i 

1 l^M,N=l,i-y5 j/-ni J 



A^?"[ M+W —xi—m ^ k' —m ^\ r/—m 



il j/—m J (2ir) d (2ir) d ty—m <-ty+k'' -tf-m If —m ' j/- 



(157) 
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